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We calculate the exchange energy in two dimensional laterally coupled quantum dots using Heitler- 
London, Hund-Mullikan and variational methods. We assess the quality of these approximations 
in zero and finite magnetic fields comparing against numerically exact results. We find that over 
the much simpler Heitler-London method surprisingly, the Hund-Mullikan method does not offer 
any significant improvement, whether at large or small interdot distances. Contrary to that, our 
variational ansatz proves substantially better in a single dot at finite magnetic field, all approximate 
methods fail. This reflects the qualitative change of the single electron ground state from non- 
degenerate (harmonic oscillator) to highly degenerate (Landau level). However, we find that the 
magnetically induced failure does not occur in the most important, double-dot, regime. 



I. INTRODUCTION 

In the last decade, a great interest in quantum dots 
has aroused, due to their potential use as hardware for 
a scalable quantum computer [SJ [TU 031 [T7] . In this im- 
plementation, the electron spin in quantum dots is used 
as the basic unit of information (qubit) [13 . To per- 
form any quantum algorithm, 1-qubit and 2-qubit gates 
are sufficient. Any 2-qubit gate can be achieved through 
single qubit rotations and an adequate switching of the 
exchange energy, which parametrizes the spin coupling in 
the Heisenberg spin exchange Hamiltonian [JJ [TT| . The 
knowledge of the exact value of the exchange energy is 
important as it determines the time of the the \J SWAP, 
a fundamental 2-qubit gate (TTJ. The properties of the 
exchange energy in lateral dots have been investigated 
by a variety of methods: Heitler-London [31 HI 151 IH1 |2"4] . 
Hund-Mulliken [H [H], Molecular Orbital [H], Varia- 
tional [16) . Configuration Interaction |20 [ l2T | [32|. Hartree 
[25) . Hartree- Fock [141 l2"5j, Hubbard model JT3J, quantum 
Monte Carlo [S3] and local spin density approximation 



This work compares several standard methods 
(Heitler-London, Molecular Orbital, Hund-Mulliken, 
Variational and the exact Configuration Interaction) to 
compute the exchange energy in quantum dots. We find 
that, with the exception of our Variational ansatz, the 
other standardly used extensions of the Heitler-London 
method do not in fact offer any real improvement. We 
revisit the failure of the Heitler-London method in a fi- 
nite magnetic field [¥|. We explain the failure as due to 
the qualitative change of the single electron ground state 
from non- to highly- degenerate. As our most important 
result, we find that contrary to the singlet dot case, the 
Heitler-London method does not suffer the failure in the 
magnetic field in the double dot regime. 

The article is organized as follows. In Sec. [II] we de- 
fine the double dot model, and present the Fock-Darwin 
states. In Sec. |III| we define the approximate methods 
that we compare in further. In Sec. |IV| we present the 



comparison, varying the magnetic field and the interdot 
distance. 



II. MODEL 

We assume to have two electrons in a harmonic electro- 
static potential with one (quantum helium) , or two sym- 
metrical (quantum hydrogen) minima llb| I25| I30[ 131] . 
We consider the electron to be two-dimensional, as appro- 
priate for electrically defined lateral semiconductor quan- 
tum dots [HlIIIlinilSB] in GaAs/AlGaAs heterostructure. 
We describe the electrons using the single band effective 
mass approximation |10j . A constant magnetic field is 
applied along the growth direction (being here along the 
z-axis) [51 12"5]. 



A. The Hamiltonian 

For the quantum hydrogen (double dot) the Hamilto- 
nian of the i-th electron is [H] 
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where the kinetic momentum P = p+eA is expressed us- 
ing the canonical momentum p = — iSV, and the vector 
potential A = — ™R x B, projected to the two dimen- 
sional plane. The magnetic field B — (0,0, B), and the 
position vector R = (r,z), r = (x,y). The positron ele- 
mentary charge is e, the effective mass of the electron is 
m, and hu is the confinement energy [31 . The vector d 
defines the main dot axis with respect to the crystallo- 
graphic axes. 

The quantum helium (single dot) can be described set- 
ting the interdot distance to zero [9], resulting in the 
Hamiltonian, 
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For a system of two electrons the Hamiltonian is 
H = H { d 1) +H { d 2) +H c , 



(3) 



where He is the Coulomb interaction Hamiltonian be- 
tween the two electrons, 



H c = 
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where eo is the vacuum dielectric constant and e r is the 
dielectric constant of the material. We neglect the Zee- 
man and spin-orbit interactions |27j , as they are substan- 
tially smaller than the above terms. 

In numerical computations, we use the parameters of 
GaAs: e r = 12.9, m — 0.067m e (m e is the free electron 
mass), take the confinement energy hui = 1 meV |28j . 
and place the dot such that d = (d, 0), that is along the 
x-crystallographic axis. 



B. Fock-Darwin states 

The eigenfunctions of the single dot Hamiltonian, 
Eq. Q , are the Fock-Darwin states |5J > 

V£?(r) = C ne r\%-^^L^ (£) e>r (5) 

Here C n i is a normalization constant, L l n are Laguerre 
polynomials, and n, I are principal and orbital quantum 
numbers, respectively. The right hand side of Eq. Q is 
expressed in polar coordinates (r, ip) . The corresponding 
energies are 
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is the effective confinement length. 



C. Exchange Energy 
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where Si,2 = fi.<T 12 /2 are spin operators and the only pa- 
rameter is J, the exchange energy. By this construction, 
the exchange energy is defined as the difference between 
the energy of the triplet and the singlet state [T3] . There- 
fore, we can write |22| : 



(9) 



where tyg and are the spin singlet and triplet two- 
electron wave functions, respectively. Below, these wave- 
functions are constructed from the single electron states 
(Fock-Darwin) appropriately, according to the specific 
method. After that, we evaluate integrals in Eq. |9]) nu- 
merically, to obtain the exchange energy. 



III. METHODS 

In this section we define four approximative methods 
to compute the exchange energy for a system of two elec- 
trons in a coupled double dot system. They differ in 
the way how the two electron states entering Eq. ([9| are 
constructed. 



A. Heitler-London method 

The Heitler-London approximation is the simplest 
method to calculate the exchange coupling in a two elec- 
tron dot [51 Q3]. It uses the lowest Fock-Darwin state 
(n = £ = 0), displaced to the potential minima positions 
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From these one constructs the two-electron singlet and 
triplet states as 

m s/T = C s/T [tp L (ri)ip R (r 2 )±^R(ri)^ L (r 2 )}xs/T, (H) 

where the singlet \s = (l/V^) (| t ® 4- - i ® t)) 5 
and xt is one of the three triplet states, 
(l/V2)(|t®4- + + ®t», | t ® t), U ® +>. and 
Cs/t is a normalization constant. 



Neglecting the Coulomb interaction, one can write the 
two-electron wavefunction using single electron eigen- 
states. If the state is separable to the spinor and or- 
bital parts, to have an antisymmetric function, orbital 
part must be symmetric and spinor antisymmetric, or 
vice versa. At low-temperature the relevant two electron 
Hilbert space can be restricted to comprise the two lowest 
orbital eigenstates of the Hamiltonian in Eq. (|3|. Such 
restricted subspace is described by an effective Hamilto- 
nian 



B. Molecular Orbital method 

The basic idea of the Molecular Orbital method is 
to use molecular orbitals, instead of the localized Fock- 
Darwin states, as the basic building blocks of the two 
electron wavefunctions |5J E] • Molecular orbital is the 
single electron eigenstate of the double dot Hamiltonian, 
Eq. 0. We take the following one-electron wavefunc- 
tions as approximations to the lowest two molecular or- 
bitals 
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where C± are the normalization constants. The sym- 
metrized combinations of the previous wavefunctions 
form the following two electron states, 



*i = hMri)V+(r 2 )ks, 

* 2 = [V>-(ri)^_(r a )]xs, 

* 3 = C 3 [^+(ri)V-(r 2 ) + V+(r 2 )^-(ri)] xs, 

* 4 = C 4 [^+(ri)V-(r 2 ) - ^+(r 2 )V-(ri)] Xt- 



(13) 



Here again, C3 and C4 are normalization constants. 
To obtain the two electron energies, we diagonalize the 
Hamiltonian in Eq. ^ in the basis consisting of these 
four states. 



C. Hund-Mullikan Method 

In this method one expands the two electron basis used 
in the Heitler-London method, Eq. ( fll] ), by doubly oc- 
cupied states @], ipL^i^L^) and ipR(ri)^ R (r 2 ). It is 
simple to check that due to the choice of the single elec- 
tron molecular orbitals, Eq. (12 1, such expanded basis is 



equivalent to the one used in the Hund-Mullikan method. 
It then follows that the Hund-Mullikan approximation is 
equivalent to the Molecular Orbital with the single elec- 
tron orbitals approximated by Eq. (12). 



D. Variational method 

In the variational method [TB] one makes a guess for 
a trial wavefunction, which depends on variational pa- 
rameters. These parameters are adjusted until the en- 
ergy of the trial wavefunction is minimal. The resulting 
trial wavefunction and its corresponding energy are vari- 
ational method approximations to the exact wavefunc- 
tion and energy. Here we use the Heitler-London type of 
ansatz, Eq. (11), with 
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The energy is defined as the minimum over the varia- 
tional parameter D, 
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where H is given in Eq. (|3]) and the minimization is done 
numerically 



E. Configuration Interaction Method (Exact) 

The configuration interaction is a numerically exact 
method |32) , in which the two electron Hamiltonian is di- 
agonalized in the basis of Slater determinants constructed 
from numerical single electron states in the double dot 
potential [TJ|2j|32]. Typically, we use 21 single electron 
states, resulting in the relative error for energies of order 

icr 5 0. 



IV. RESULTS 

Here we present results of the numerical calculations 
for the exchange energy using the methods listed in the 
previous. In Fig. [I] the exchange energy is plotted as 
a function of the interdot distance for two electrons in 
a double dot in zero magnetic field. One can see that 
for large interdot distance the exchange energy falls off 
exponentially, a fact that all methods reflect correctly. 
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FIG. 1: Exchange energy as a function of the interdot distance 
for two electrons in a double dot in a zero magnetic field. 

This suggests that, at least in principle, an efficient 
control of the exchange energy can be achieved by in- 
creasing the potential barrier separating the dots and/or 
by increasing the interdot separation. At a small inter- 
dot distance all methods differ significantly from the ex- 
act result. We note here that, despite the common be- 
lief, the Hund-Mullikan (equivalent to Molecular-Orbital) 
method does not offer any improvement, even at small 
interdot distances (strong interdot couplings) . The Vari- 
ational method of the form that we choose, on the other 
hand, proves more robust, typically cutting the error of 
the Heitler-London to a half. We will see on examples 
that follow, that these two features are generic and we 
explain them below. Figure [2] shows the exchange energy 
as a function of the interdot distance in a finite magnetic 
field IT. The exchange energy decreases faster than in 
zero magnetic field. A simple explanation follows from 
noting that the natural length scale is the effective con- 
finement length £b, which drops with the magnetic field, 
as seen from Eq. 0. 

Figure [3] shows the exchange energy as a function of 
the magnetic field for a single dot. We can see that as 
the magnetic field increases, the approximative methods 
become increasingly off the exact results. Worse than 
that, except the variational method, even the trend is 
wrong (growth, instead of a fall-off) 

To understand this failure, we plot in Fig. [4] the sin- 
gle electron single dot (Fock-Darwin) spectrum. One 
can appreciate the qualitative change between the low 
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FIG. 2: Exchange energy as a function of the interdot distance 
for two electrons in a double dot in a finite magnetic field of 
1 Tesla. 
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FIG. 3: Exchange energy as a function of the magnetic field 
for two electrons in a single dot. 
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FIG. 4: Spectrum of a single electron in a single dot (Fock 
Darwin) as a function of the magnetic field. 



and high field limit. In the first, the ground state is 
non-degenerate, while in the second, a highly degenerate 
Landau level forms. The higher degeneracy, the more 
mixing and therefore worse results are expected for all 
methods based on a basis built from just a few single 
electron states. In other words, the effective strength of 
the Coulomb interaction grows with diminishing of the 
energy separation of the single electron states due to the 
magnetic field |151 ITS] . 

To gain further insight, we plot in Fig. [5] a comparison 
of the singlet and triplet energies in the Heitler-London 
approximation with their exact counterparts as a func- 
tion of the magnetic field for a single dot. From this we 
can see that both the singlet and the triplet are in the 
Heitler-London approximation similarly off the exact re- 
sults (that is, the failure is not due to just one of them). 
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FIG. 5: Singlet and triplet energies in the Heitler-London 
method with their exact counterparts for a single dot. 
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FIG. 6: Exchange energy as a function of the magnetic field 
for two electrons in a two-electron double dot (d = 60 nm). 

If operating as a qubit, the double dot will be ma- 
nipulated at large interdot distances. Figure [6] shows 
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the error of the singled is reduced, while the triplet is 
untouched. Since the exchange energy is the difference 
of the two energies, using the Hund-Mullikan/Molecular- 
Orbital is in fact detrimental compared to the Heitler- 
London method. The Variational method that we choose, 
on the other hand, treats the symmetric and antisymmet- 
ric wavefunctions similarly. It is then natural to expect 
that their errors tend to compensate, resulting in more 
precise value for the exchange energy. 
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FIG. 7: Singlet and triplet energies in the Heitler-London 
method with their exact counterparts for double dot (d — 60 
nm). 



CONCLUSIONS 



the exchange energy in this regime. Surprisingly, the ap- 
proximate methods reflect the exponential supression of 
the exchange energy with the magnetic field correctly. 
This can be traced down to the fact that in the Heitler- 
London method, the exchange energy is proportional to 
the overlap of the localized single electron states, from 
Eq. (10 1. This is illustrated further in Fig. [7] where we 
plot the singlet and triplet energies in the Heitler-London 
approximation together with their exact values in this 
regime. One can see that even though the approximate 
values do not converge to the exact one, their errors tend 
to compensate. As a result, the exchange energy falls 
towards zero as the magnetic field is enlarged. This find- 
ing also helps to understand the previously seen trends. 



Namely, as seen from Eq. (13 1, the Molecular-Orbital ap- 
proximation expands, compared to the Heitler-London, 
the singlet subspace only. In the resulting energy pair, 



To conclude, we studied the exchange energy in two 
electron single and double lateral quantum dots. We 
compared four approximate methods: Heitler-London, 
Hund-Mullikan, Molecular Orbital and Variational, with 
numerically exact results (the configuration interaction 
method). We find that, compared to the much simplet 
Heitler-London method, the Hund-Mullikan and Molec- 
ular Orbital methods do not offer any improvement, 
whether at large or small interdot distances. We explain 
that noting the former two methods treat the singlet and 
triplet wavefunctions differently leading to uncompen- 
sated errors. On the other hand the variational ansatz 
proves robust. At finite magnetic field, all approximate 
methods fail badly. This is a consequence of the qualita- 
tive change of the single electron ground state. Finally, 
and most important, we find that all the approximate 
methods we study are free from the failure in the double 
dot regime. 
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